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ABSTRACT 



Context. Particle acceleration consequences from fluctuating electric fields superposed on an X-type magnetic field in 
coUisionless solar plasma are studied. Such a system is chosen to mimic generic features of dynamic reconnection, or 
the reconnective dissipation of a linear disturbance. 

Aims. Explore numerically the consequences for charged particle distributions of fluctuating electric fields superposed 
on an X-type magnetic field. 

Methods. Particle distributions are obtained by numerically integrating individual charged particle orbits when a time 
varying electric field is superimposed on a static X-type neutral point. This configuration represents the effects of 
the passage of a generic MHD disturbance through such a system. Different frequencies of the electric field are used, 
representing different possible types of wave. The electric field reduces with increasing distance from the X-type neutral 
point as in linear dynamic magnetic reconnection. 

Results. The resulting particle distributions have properties that depend on the amplitude and frequency of the electric 
field. In many cases a bimodal form is found. Depending on the timescale for variation of the electric field, electrons 
and ions may be accelerated to different degrees and often have energy distributions of different forms. Protons are 
accelerated to 7-ray producing energies and electrons to and above hard X-ray producing energies in timescales of 1 
second. The acceleration mechanism is possibly important for solar flares and solar noise storms but is also applicable 
to all coUisionless plasmas. 
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1. Introduction 

, Release of stored magnetic energy via particle acceleration 
■ is a characteristic feature of astrophysical plasmas. In the 
[ particular case of the Sun, we see this manifested in the 
. catastrophic events of flares, as well as in quieter phenom- 
' ena like radio noise storms. Similar phenomena are observed 
in other late- type stars, and similar physics may be involved 
I in understanding a wide variety of astrophysical objects (see 
' e.g. Kuiipers 119931 Hanasz and Lesch [20051) . 

The special case of solar flares involves particular chal- 
lenges to theory. A large fraction (several tens of percent) 
of the flare energy is manifested initially in the form of 
fast electrons (accelerated out of the background distribu- 
tion to '--^lOO keV in about 1 second and to ^ 100 MeV 
in a few seconds), which reveal their presence by produc- 
ing bremsstrahlung X-rays (e.g. Miller 1998^ MacKinnon 
I2006p . Protons are accelerated in flares to energies of sev- 
eral tens of MeVs in a timescale of one second (Miller 
119981 Aschwanden l2002| ). Thus the acceleration of particles 
is an important part of the energy release process, rather 
than an energetically unimportant consequence of the flare. 
Moreover, radio signatures (Type I noise storms, Type III 
bursts away from flares) testify to particle acceleration at 
"quiet" times. 

Magnetic reconnection is one of the primary candidate 
mechanisms for releasing non-potential energy from mag- 
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netized plasmas (e.g. Priest and Forbes '200 f). The electric 
field in the current-carrying region also makes it a natural 
particle accelerator. Collision-dominated sheets will involve 
the production of some runaway particles, but almost by 
definition particle acceleration is not a primary consequence 
of such a situation (e.g. Smith I1980p . However, Martens 
()1988p gave order-of-magnitude arguments in favor of a col- 
lisionless current sheet as both the energy release mecha- 
nism and the particle accelerator in flares. Particle accel- 
eration is energetically the primary consequence of such a 
situation. CoUisionless reconnection thus assumes great po- 
tential importance in understanding the flare process, par- 
ticle acceleration, energy conversion and release in astro- 
physical plasmas generally (Petkaki and MacKinnon I1997| 
Heerikhuisen et al. i2002i Hamilton et al. 2003; Turkmani et 
al. 2006; Wood an d Neu kirch [2005 Vainchtein et al. [2005} 
McClements et al. [2006]) . 

Here we present test particle calculations designed to 
illuminate the consequences for particle acceleration of dy- 
namic reconnection. We have in mind particularly the pic- 
ture of Craig and McClymont (1991, 1993), in which a lin- 
ear disturbance passes through a magnetic configuration 
containing an X-type neutral point. The disturbance trav- 
els non-dissipatively with the local Alfven speed until it 
approaches the dissipation region surrounding the neutral 
point, where the resistive diffusion term in the induction 
equation becomes important. The wave damps resistively 
in a few system transit times, with consequences (heating 
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or particle acceleration) determined by the physical nature 
of the resistivity. Several assumptions were made in Craig 
and McClymont's original discussion (linear disturbance, 
cold plasma limit, 2-D, Ohm's law including only a scalar 
resistivity) but this essential picture still affords a qualita- 
tive guide in more complex situations (e.g. McClymont and 
Craig fT996l Senanayake and Craig l2006|) . 

In Petkaki and MacKinnon (|1997p . we examined the be- 
havior of protons in the presence of electric and magnetic 
fields obtained from the Craig and McClymont ()1991l) anal- 
ysis. Here we carry out a complementary exercise, studying 
test particle evolution in the presence of simple fields chosen 
to mimic generic features of dynamic reconnection. Our aim 
is to comment on particle acceleration consequences, in a 
parametric way that does not depend on a particular set of 
simplifying physical assumptions or boundary conditions. 
Time-dependence of the electric field is the essential ingre- 
dient reflecting the dynamic character of the reconnection. 
We present examples of distributions resulting from a time- 
independent electric field for comparison and highlight dis- 
tinct features of the distributions resulting from dynamic 
situations. 

We use the Craig and McClymont (|199ip linear solution 
as a qualitative guide for the spatial and temporal form of 
the electric field. Our adopted field also resembles a linear 
situation in displaying a time dependence that does not 
change (i.e. does not develop multiple frequencies, saturate, 
etc.). Basing our calculations on this linear picture makes 
it unlikely that they will provide a complete description of 
what happens in a flare, although they offer useful insight. 
They may however be particularly relevant to non-flaring 
particle acceleration, e.g. in solar noise storms, or as part 
of the explanation of 'quiescent' radio emission seen from 
RS CVn binaries (Kuijpers & van der Hulst fTMS)) . 

Since we aim to emulate a linear situation we may pick 
our test particles from an isotropic, homogeneous distribu- 
tion representing the background. This is in contrast to par- 
ticle studies of nonlinear reconnection, where consistency 
demands consideration of the motion of particles into the 
dissipation region. In most studies particles are injected in 
two opposite quadrupoles of the X-point and they subse- 
quently are driven, by the Ex B drift due to an imposed con- 
stant electric field, to cross the nonadiabatic region or miss 
it depending on their initial conditions (see e.g. Burkhart 
et al.[Il90). 

Many previous studies of test particle evolution in 
steady reconnection exist. Here we mention particularly the 
work of Martin (1986), which demonstrates that the orbits 
of such test particles are chaotic, and of Burkhart et al. 
(|19901 I199ip who iterated from the test particle calcula- 
tions to construct a self-consistent description of the dif- 
fusion region. Recent work studies regular and chaotic dy- 
namics in 3-D reconnecting current sheets (Efthymiopoulos 
et al. 120051 Gontikakis et al. I2006P or studies particle or- 
bits in the presence of 3-D magnetic nulls (Heerikhuisen 
et al. ?IUU^ Dalla and Browning [MJ5)) . Particularly rele- 
vant here is the exploratory, analytical study of Litvinenko 
(|2003p which looks at charged particle orbits in an oscil- 
lating electric field in a magnetic field containing a neutral 
line. 

The next section gives details of the specific model 
we adopt in order to study particle acceleration in time- 
dependent reconnection, while Section [3] describes our re- 



sults for particle distributions. Section[4]discusses some pos- 
sible implications of our results. 

2. Model for Particle Acceleration in fluctuating 
electric fields 

We are going to study the evolution of test particles in the 
presence of electromagnetic fields chosen to mimic generic 
features of dynamic reconnection. Time-dependence of the 
electric field reflects the dynamic character of the reconnec- 
tion. 



2.1. Equations of motion 

We solve numerically the relativistic equations of motions of 
test particles (particles are expected to acquire relativistic 
velocities) in electromagnetic fields and in the observer's 
reference frame: 



dt 



P 

7717 



dp 1 
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(1) 
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where 7 = (1 — (it/c)^)^"'"/^, u = ^p/rwy. 

To model the reconnection magnetic field, we adopt an 
idealized 2-D magnetic field containing an X-type neutral 
point: 



B 



So 
D 



(yx + xy). 



(3) 



The current density vanishes for this field configuration. 
The field lines are the solutions of || = ^ which are hyper- 
bolae y"^ — x"^ — const. The X-line (neutral line) lies along 
the z-axis. The field strength depends on position thus: 



(4) 



where x'^ + y"^ — . Note that this configuration has no 
natural scale length. Requiring the field to have a value of 
10^ gauss at a typical active region distance of 10^ cm from 
the neutral point, fixes only Bq/D = 10~^ gauss cm^^. We 
are free to use other considerations to fix one of Bq and D 
independently, as we do below in introducing dimensionless 
variables. An electric field is imposed in the z direction, 
with spatial and temporal form chosen to mimic qualitative 
features of dynamic reconnection (see Sec. 12. 2p . 

We normalize distances to Dn and times to the gyrope- 
riod at r — Z?„. We denote the resulting timescales by Tp 
and Te for the cases of electrons and protons respectively. 
As noted above, D„ is as yet undetermined. It turns out to 
be convenient in this relativistic calculation to choose Dn 
such that velocities are normalized to the speed of light. 
This has the consequence that £>„ takes different values 
De = Cy/{nieD/eBo) and Dp ~ Cy/{mpD/eBo) for elec- 
trons and protons respectively (Petkaki and MacKinnon, 
1994, Petkaki. such that 



„ . 77l„ . 1 „ 

Dp = {^)-D, 

777e 



(5) 



Specifically, with Bo/D = 10"'^, we find D^ = 1.3 x 10^ cm 
and Dp = 5.6 x 10^ cm. With our choices of E and B the 
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Lorentz equations ([T]) and ^ become in three dimensions 
and in dimensionless units: 



dx 

Is 

dy 
dt 
dz 

dpx 
di 

dpy 
dt 

dpz 
dt 



7 
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^yuz 

E + £{xux + yuy) 



(6) 



where 7 = (1 +^1^ ^ Py^ +Pz'^)^ and e = +1 for protons, 
e = — 1 for electrons. E{r,t) is the true value of the elec- 
tric field and E = E{D / Bo)Di is the dimensionless electric 
field, with subscript i taking the values e for electrons or p 
for protons. Energies are now normalised to the particle rest 
mass energy so that kinetic energy in dimensionless units 
is just Kkin =7 — 1. Equations ^ with appropriate initial 
conditions and a specific form for E describe the motion of 
a particle. 

2.2. Electric Field 



Craig and McClymont (|199ip guide us in adopting a func- 
tional form of electric field which allows us to investi- 
gate consequences of time-dependence in a parametric way. 
Their resistively damping, linear disturbance involves a reg- 
ularly oscillating electric field whose amplitude is greatest 
in the region near the neutral point where the resistive term 
of the induction equation becomes important. Far from this 
region the disturbance is Alfv'enic in character and domi- 
nated by the boundary conditions, so that the electric field 
amplitude always maximises in the central, diffusion region 
(see also Petkaki and MacKinnon 119971 Section 3). Thus 
we adopt the following form for the electric field E: 



E EQsa\{ujt)zf{x,y) 



(7) 



where f{x,y) describes the spatial variation of E- We take 



fix,y) = exp{-a^y/{\r\)) 
where, ar, — 2.5 x 10^^, a„ 



(8) 



3.776x10-2, 



-y 



Consistent with our concentration on particle acceleration 
near the neutral point and with in situ measurements 
in Earth's magnetosphere (0ieroset et al., I200ip . we ex- 
pect that resistivity will be primarily inertial in character 
(Speiser I1970p . As shown in Figure [1] this form approxi- 
mates the radial form of the electric field calculated from 
the Craig and McClymont (1991) solution, for an (iner- 
tial) resistivity 77 estimated assuming the dominant contri- 
bution from 1 keV protons (see Petkaki and MacKinnon 
fT997[ Section 3; Speiser, \T9Wf . The form of the electric 
field is shown in Figure [T] at t = 0. The exact Craig and 
McClymont solution develops more complex spatial struc- 
ture but the spatial form of Eq. ([8]) thus embodies a dissi- 
pation region on the appropriate length scale. 




1.00 

r (Dp) 



Fig. 1. Electric field (solid line) calculated from the Craig 
and McClymont p99ip solution for inertial resistivity rj — 
3.1724 X IQ-^^ and for the fimdamental mode of azimuthal 
symmetry (n=0) at t=0. Approximate electric field de- 
scribed in Equation [S] is shown in dotted line at t=0. 

The frequency of oscillation of the electric field which 
we denote by w is a free parameter. Each simulation uses 
one value of w. We take values of uj such that 1/1000 < 
to < 10000, corresponding to a broadband wave spectrum 
which may propagate in such a system (see Petkaki and 
Mackinnon [TM7|) . 

To compare with a simple, unvarying state, we also cal- 
culate the energy distributions that result from a constant 
in time imposed electric field 

^^Eozf{x,y) (9) 

where f{x,y) is defined in Equation [H] With its nonzero 
curl, this assumed form of E cannot represent a steady 
state reconnection. We employ it primarily to provide a 
simple, unvarying state for comparison with results in the 
time-dependent situation. 

Key to understanding particle behaviour near the neu- 
tral point is the 'adiabaticity' radius Vad, the distance from 
the neutral point at which the Larmor radius equals the 
magnetic field scale length. For r > rad, particles move adi- 
abatically. In other words, if the distance of the particle 
from the neutral point is of the order of its Larmor radius, 
then the particle is non-adiabatic. The 'adiabaticity' radius 
depends on the particle mass and velocity perpendicular to 
the magnetic field, u±, and is given by 

(I^^)i (10) 



rad 



eBo 



For electrons and protons of the same energy, the electron 
gyroradius is (^)^ smaller than the proton gyroradius. 

2.3. Numerical Method 

Due to the complexity of the orbits, their calculation cannot 
be done analytically. For integrating the ordinary differen- 
tial equations (ODEs) describing the motion of the parti- 
cles, we use the Bulirsch-Stoer method (Press et al. I1996p . 
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For a single particle with the same initial conditions, the 
orbit changes if the accuracy required of the integration 
routine is varied, and when the particle crosses the neutral 
point area (for general properties of X-type neutral point 
orbits see e.g. Martin fTOSGp . The statistical properties of the 
distribution of test particles, which are of primary interest 
here, are unaffected by changes in the accuracy required of 
the integration routine. In the absence of an electric field, 
the routine conserves particle energy to one part in 10~^. 

We start the integration of particle orbits at t=0 and 
with the particles positioned randomly in a box with the 
following size 



1.0 < xO < 1.0 
-1.0 < yO < 1.0 
zO = 0.0 



(11) 
(12) 

(13) 



in electron or proton units depending on the species. We 
integrate the particle orbits up to 230400 timesteps (rg) for 
electrons and 5360 (rp) for protons. With Bq/D = 10~^ and 
our form of dimensionless units these times correspond to 
1 second for electrons and protons. The initial velocities of 
the particles are picked randomly from a Maxwellian dis- 
tribution of temperature 5 x IQ^K (~ 431 eV),_a typical 
coronal value. We consider only small values for Eq, consis- 
tent with the passage of a disturbance in the linear regime 
(Craig and McClymont, 1991). Values of 0.0001, 0.001 and 
0.01 are used in the actual calculation. The value 0.001 cor- 
responds to electric field = 5.88 x 10~^ statvolt/cm. These 
are moderate values for electric fields present in the solar 
atmosphere (see Foukal et al.. il986|) . 

2.4. Particle Orbits 

We are going to examine a typical proton orbit which is 
shown in Figure [2] The amplitude of the time- varying elec- 
tric field is Eq = 0.001 and the frequency is w = 0.2. In 
Fig. [2^ we plot the (dimensionless) energy of the proton 
as a function of time for the interval 2400-5360 Tp. In Fig. 
[^b is shown the projection of the same orbit on the X-Y 
plane and in Fig.[2t the projection of the same orbit on the 
X-Z plane. Our model assumes a system scale of ~ 178Dp, 
corresponding to a typical active region scale of 10^ cm. 
In Fig. we zoom in close to the neutral point to ob- 
serve the particle orbit in detail and we look in an area of 
-l.ODp < X < l.ODp and -l.ODp < y < l.ODp. In the 
same figure superposed in dotted line are some of the mag- 
netic field lines showing the structure of the X-type mag- 
netic neutral point. All field lines tend to the separatrices 
(shown in dashed lines) as the distance from the neutral 
point becomes very large. 

The thick solid circle has radius d ~ IQrad for a thermal 
proton. Inside this region the gyroradius (Larmor radius) 
of most particles is not well defined since the particle is not 
bound to one magnetic field line and meandering motion 
is observed. The electric field accelerates or decelerates the 
proton causing further changes in the particle gyroradius 
and energy. This behavior resembles a stochastic-type ac- 
celeration. Stochasticity is introduced by the phase of the 
electric field and the phase of the particle orbit and is sus- 
tained because of the form of the magnetic field (e.g. Martin 
I1986p . Outside the magnetic neutral point area the particle 
is moving along a particular magnetic field line. The gyrora- 
dius and the parallel velocity decrease as the particle moves 



away from the neutral point. The particle mirrors and re- 
crosses the non-adiabatic region and the process is repeated 
until the end of the integration time or until the particle 
escapes the outer boundary of the system {x,y > 178Dp). 

Looking back at Fig. [5^ we see the variation of the par- 
ticle energy as a function of time. The intervals of energy 
conservation correspond to the times the particle is away 
from the nonadiabatic region, and in regions where the elec- 
tric field is decreasing. Changes in the particle energy take 
place during the crossing of the non-adiabatic region. 

The orbits of electrons show similar properties to that 
of the protons. Additionally from the way we pick our ini- 
tial conditions the electrons start closer to the neutral point 
(see Sect. 12. 2p . But since the form of the electric field for 
both species is calculated assuming proton inertial resistiv- 
ity, electrons see an area much larger than their adiabatic- 
ity radius where the electric field is close to its maximum 
value. Consequently some of the electrons start their motion 
outside their adiabaticity radius and their motion is imme- 
diately adiabatic. In this case the presence of the electric 
field does not increase their energy except if they eventually 
cross the non-adiabatic region. 

The amount of acceleration that particles get depends 
on the time they spend close to the neutral point, on the 
phase of the orbit and on the frequency of the electric field. 
We define crossing time as the time the particle needs to 
cross the non-adiabatic region (Sec. 12. 2p and is given to 
order of magnitude by 



tr. 



(14) 



where Ux,y is the velocity projection in the x-y plane. So, 



t. 



,2.828cD''^m 



1/2 3/4 



' eBo ' E^l^' 
It turns out that particles with the same energy satisfy 



(15) 



tcr^ =tcrA — )^ = 280te 



(16) 



where tcrp is the proton crossing time and tcr^ is the elec- 
tron crossing time. One would expect that in order to get 
particles effectively accelerated (or decelerated since the 
sign of the electric field is not constant) we need tE ^ tcr 
where ts = 1/w is the period of fluctuation of the electric 
field. Thus we see potential for differences between electron 
and ion acceleration. 



3. Energy Distributions of Accelerated Electrons 
and Protons 

We calculate the kinetic energy of each particle up to max- 
imum of a 1 second real time along with final positions and 
velocities. Particles in our calculation spend a relatively 
short time close to the neutral point but they get trapped 
in the magnetic configuration and re-cross the neutral point 
a number of times (Sect. [^^ . Particles encounter the non- 
adiabatic region, the process resulting in a Fermi-type ac- 
celeration. A similar phenomenon has been noted for a mul- 
tiple neutral point configuration by Kliem (1994) , in the be- 
havior of protons in the presence of an MHD disturbance by 
Petkaki and Mackinnon (|1997p . and was explored analyti- 
cally for time- varying electric field by Litvinenko (|2003p . 
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Distributions at 5360.0 Tp = 1 s 
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Fig. 3. Proton distributions for different frequencies of the electric field. The magnitude of the electric field is Eq — 0.001. 
The total integration time is 5360. 



We sample the electric field frequency range 0.001 < 
uj < 100 for protons and 0.005 < uj < 1000 for electrons. 
In Fig. [3] we plot histograms of the logarithm of the ini- 
tial and final energy distributions of protons, for magni- 
tude of the electric field Eq = 0.001 and total number of 
timesteps 5360. Each distribution is generated using 50000 
test protons. We also calculated the distributions resulting 
when the magnitude of the electric field is Eq ~ 0.01 and 
Eq — 0.0001. The initial Maxwellian distribution is shown 



in dotted lines in each panel. Panel (1) shows the distri- 
butions for electric field magnitude Eq = 0. We observe 
no change in the form of the distribution since no acceler- 
ation is taking place (see also Sect. 12. 3|) . Panel (2) shows 
the energy distributions for constant electric field (w = 0) 
to provide a comparison with steady state magnetic recon- 
nection. The final energy distribution for lu ^ has two 
distinct peaks, one at the initial Maxwellian distribution 
and a beamlike distribution close to Kkin — 10 MeV. The 
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energy distributions for constant electric field are not power 
laws as found elsewhere (see e.g. Bulanov and Sasarov ri976l 
Bruhwiler and Zweibel 1992). In those former calculations 
the particles crossed the neutral point only once, whereas 
particles recross the neutral point numerous times in our 
model. The frequency of the time-varying electric field in- 
creases progressively from cu — 0.001 in panel (3) to w = 10 
in panel (12). For w = 0.001 a small beamlike structure 
appears at Kkin = 0.3 MeV. The final energy distributions 
are bi-modal from uj = 0.01 to w = 10. For higher frequen- 
cies the proton energy distributions do not show significant 
energy changes. 

In Fig. m we plot the mean of the logarithm of the initial 
and final proton energy distributions versus the frequency 
of the electric field and for three amplitudes of the electric 
field Eo = 0_.0001 (dashed star line), Ea = 0.001 (solid star 
line), and Eq = 0.01 (dotted star line). The mean energy 
for the constant electric field is represented on this plot by 
UJ — 10^^. The same representation is used in Fig. [5] where 
we plot the standard deviation of the logarithm of the initial 
and final proton energy distributions versus lu for the same 
three amplitudes of the electric field as in Fig. IH We use 
the mean value of the logarithm of the energy to better 
represent the changes in highly non-thermal distributions. 

The mean energy increases monotonically with Eq for 
the constant electric field case and for all frequencies of the 
electric field except for the highest frequency used in our 
model. The highest mean energy is achieved for constant 
electric field for all values of the electric field. For Eq = 
0.0001 the highest energy gain for the time- varying electric 
field is achieved when 0.2 < u; < 2.0, indicating a resonant 
acceleration process. For Eq = 0.001 a peak in the mean 
energy is also present when 0.2 < oj < 2.0. For Eq — 0.001 
protons gain most energy from the low lu electric field (see 
Fig. 2]). When Eq = 0.01 protons a peak in the mean energy 
is present when 0.01 < u) < 2.0. For frequency w = 10 and 
greater the energy distribution does not change significantly 
for all values of Eq. 

In Fig. [S] we plot histograms of the logarithm of the 
initial and final energy distributions of electrons, for mag- 
nitude of the electric field Eq = 0.001 and total num- 
ber of timesteps 230480. Again we generate each distri- 
bution using 50000 test electrons, and show the initial 
Maxwellian distribution in dotted lines in each panel. Panel 
(1) shows the energy distributions for non- varying electric 
field (w = 0). The final energy distribution (shown in solid 
line) includes a small beamlike component. The lower en- 
ergy part of the final distribution is Maxwellian-like peaking 
at kinetic energy Kkin = 10~^ MeV, with a small beamlike 
component superposed at around kinetic energy Kkin — 1 
MeV. The frequency of the time-varying electric field in- 
creases progressively from uj = 0.001 (panel 2) to w = 500 
(panel 12). The bulk of the distribution is accelerated for 
the frequency range uj = 0.001 to uj — 50. Accelerated dis- 
tribution for OJ = 0.001 (panel 2) is Maxwellian-like peak- 
ing around ~ Kkin = 10~^ MeV. From uj = 0.01 (panel 
3) to w = 1 (panel 6) the accelerated distributions have 
maximum close to ~ Kkin = 0.01 MeV with maximum 
energy close to 1 MeV. From uj = 5 (panel 7) to w = 20 
(panel 9) the accelerated distributions have maximum close 
to ^ Kkin =0.1 MeV with maximum energy exceeding 
to 1 MeV. The final energy distributions are bi-modal for 
ui = 50 to u! = 500 containing a Maxwellian-like part at the 
energy range of the initial Maxwellian distribution and an 



Table 1. Percentage of Accelerated Protons for Eq = 0.001 



UJ 


> 1 


> 10 


> 100 


> IMeV 


(MeV) 


0.0 


89.1 % 


50.7 % 


14.5 % 


3.0 % 


15.9 


0.001 


59.7 % 


12.5 % 


2.9 % 


0.01 % 


1.31 


0.01 


47.9 % 


24.3 % 


8.8 % 


0.9 % 


5.86 


0.05 


51.0 % 


19.0 % 


11.1 % 


2.9 % 


8.0 


0.1 


55.2 % 


22.9 % 


14.3 % 


3.9 % 


5.9 


0.2 


52.9 % 


27.7 % 


17.8 % 


5.2 % 


5.0 


0.5 


72.7 % 


52.0 % 


37.4 % 


17.9 % 


18.5 


1.0 


82.3 % 


51.4 % 


32.6 % 


11.1 % 


9.33 


2.0 


30.9 % 


4.2 % 


2.4 % 


0.8 % 


7.33 


5.0 


6.0 % 


3.3 % 


2.3 % 


0.2 % 


3.8 


10.0 


0.3 % 


0.1 % 


0.01 % 


0.0 % 


0.31 


50.0 


0.0 % 


0.0 % 


0.0 % 


0.0 % 


0.005 


100.0 


0.0 % 


0.0 % 


0.0 % 


0.0 % 


0.005 



accelerated part at higher energies with peak in the range 
^ Kkin = 0.03 to - Kktn = 0.7 MeV. For higher frequen- 
cies the energy distributions do not show significant energy 
changes in the time of maximum 1 s. 

The electron distribution gains energy for most of the 
frequencies of the electric field that we used in this model 
except for the highest frequencies. In Fig. [7] we plot the 
mean of the logarithm of the electron energy distributions 
versus uj for two amplitudes of the electric field Eq = 0.0001 
(dashed star line), Eq = 0.001 (solid star line). In the same 
plot we superposed as error the standard deviation in the 
mean for each distribution. The mean energy for the con- 
stant electric field is represented on this plot by uj = 10~^. 
The highest energy gain achieved is for electric field with 
ti> = 50 for amplitude of the electric field Eq = 0.001. When 
Eq = 0.0001 electrons gain most energy for the constant 
electric field and for uj = 50. 

In the Tables that follow we summarized the energy gain 
aspects of the acceleration mechanism. For each particle we 
find the energy gain using its initial and final kinetic energy, 
that is: 

^K ^ final '^initial (A'7\ 

initial initial 

where Kl^^^^^^ is the initial energy of the i-th parti- 
cle and K^finai ^^'^ ?iT\Si\ energy of the i-th particle. 
AK' / Kl„^^^^i = A, where A takes the values 1, 10, 100. 

In Table [T] we list the accelerated percentage of the final 
proton distribution for each electric field frequency and for 
Eq = 0.001. The first column lists the frequencies of the 
electric field. The second, third and fourth columns list the 
percentage of the final proton distribution for which A is 
greater than 1, 10 and 100 respectively (Eq.[T7I). The fifth 
column lists the percentage of the final proton distribution 
that has energy greater than 1 MeV and the last column 
list the highest energy in the final proton distribution in 
MeV. We see that depending on the frequency of the electric 
field, ~ 0.2% to ~ 17.9% of the proton distributions get 
accelerated to 7-ray producing energies in 1 s. 

In Tables [5] and [3] we list the accelerated percentage of 
the final proton distribution for each available electric field 
frequency for Eq — 0.0001 and for Eq — 0.01. The layout 
of these tables is the same as for Table [T] discussed before. 

In Tables S] and [5] we list the accelerated percentage of 
the final electron distribution for each available electric field 
frequency for Eq = 0.001 and for Eq = 0.0001 respectively. 
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Fig. 6. Electron distributions for different frequencies of the electric_ field. The initial Maxwellian distribution is shown 
in dotted line on each panel. The magnitude of the electric field is Eq = 0.001. The maximum integration time of each 
particle trajectory is 230480. 



The first, second, third and fourth columns are layed out 
as in Table [1] The fifth column lists the percentage of the 
final electron distribution with energy greater than 20 keV 
and the last column lists the highest energy in the final 
electron distribution in MeV. For = 0.001 and for most 
frequencies of the electric field (and for constant electric 
field) the bulk of the electron distributions get accelerated 
to X-ray producing energies in the timescale of our model. 
When Eq = 0.0001 only small percentage of the electron 



distribution accelerates to X-ray producing energies except 
for constant electric field where ~ 50% accelerates to X-ray 
producing energies and for a; = 50 where ~ 23% accelerates 
to X-ray producing energies. 

Flare fast electrons as revealed by hard X-ray observa- 
tions generally have energy distributions characterised by 
energy spectral indices in the range 2-5 (Dennis, l1985|) . 
RHESSI data allow a less crude characterisation of the 
energy distribution (Kontar et al. 12005^ . but this range 
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Table 2. Percentage of Accelerated Protons for Eq = 
0.0001 





> 1 


> 10 


> 100 


(Me V ) 


0.0 


26.1 % 


4.0 % 


1.0 % 


0.36 


0.001 


9.6 % 


1.6 % 


0.2 % 


0.12 


0.01 


6.1 % 


2.3 % 


0.4 % 


0.13 


0.05 


0.14 % 


2.9 % 


0.7 % 


0.14 


0.1 


6.8 % 


4.3 % 


1.4 % 


0.17 


0.2 


12.6 % 


9.5 % 


4.6 % 


0.5 


0.5 


33.8 % 


24.0 % 


8.0 % 


0.96 


1.0 


51.1 % 


38.1 % 


15.0 % 


1.32 


2.0 


17.2 % 


11.5 % 


3.5 % 


0.21 


5.0 


2.9 % 


0.2 % 


0.0 % 


0.04 


10.0 


0.01 % 


0.0 % 


0.0 % 


0.012 


50.0 


0.0 % 


0.0 % 


0.0 % 


0.005 


100.0 


0.0 % 


0.0 % 


0.0 % 


0.005 



Table 3. Percentage of Accelerated Protons for Eq ~ 0.01 





> 1 


> 10 


> 100 


> IMeV 


(MeV) 


0.0 


100 % 


100 % 


99.2 % 


45.3 % 


525 


0.001 


92.9 % 


57.1 % 


12.8 % 


2.6 % 


24.5 


0.01 


98.4 % 


80.1 % 


38.9 % 


9.5 % 


237 


0.1 


99.8 % 


97.4 % 


67.8 % 


38.8 % 


225 


0.2 


99.8 % 


96.6 % 


69.3 % 


47.2 % 


197 


0.5 


99.8 % 


97.9 % 


84.9 % 


59.8 % 


202 


1.0 


99.6 % 


97.9 % 


88.5 % 


57.6 % 


343 


2.0 


97.1 % 


78.0 % 


15.6 % 


0.1 % 


56.5 


10.0 


40.6 % 


1.9 % 


0.04 % 


0.0 % 


0.18 


50.0 


3.0 % 


0.03 % 


0.0 % 


0.0 % 


0.005 


100.0 


0.4 % 


0.0 % 


0.0 % 


0.0 % 


0.005 



nonetheless gives a reasonable starting point for compar- 
ison with our results. With their high-energy peaks, some- 
times bimodal in form, many of the distributions shown in 
Figures[3]and [6] arc clearly some way from those implied by 
observations. The electron distributions of Figure [5] would 
give very hard photon spectra, harder than usually observed 
in flares. When segments of the distributions appeared to 
decline in energy in power-law form E^^ , we fit power-laws 
in energy to them, finding values of 6 between about 1 and 
2.5. Although we can account for particle acceleration in 
this way, to very high energies, we have to appeal to some 
other agent to redistribute energy among the accelerated 
particles to be compatible with observations. Very hard en- 
ergy distributions are also found in most other calculations 
of acceleration in either one or many dissipation regions 
(e.g. Turkmani et aL. I^IMj) . 

By fixing the duration of the integrations at 1 second 
we generate a snapshot of the distributions produced as 
particle acceleration proceeds. Obviously, particle energies 
will be less for shorter periods and greater for longer ones. 
As an illustration, in Figure [5] we plot the time evolution 
of the logarithm of the mean energy of proton distributions 
for Eq — 0.0001, as functions of ui. Mean energy is plotted 
for 0.25 s (1340 Tp), 1 s, and 2 s (10720 Tp). We observe that 
the greatest changes in mean proton energy take place for 
the lowest frequencies uj. At high frequencies, on the other 
hand, proton mean energy apparently changes little after 1 
s. A sort of steady state is approached. Since no particles 
escape, this indicates a decrease with proton energy of the 
energy increment experienced on each return to the dissipa- 



Table 4. Percentage of Accelerated Electrons for Eq = 
0.001 



LU 


> 1 


> 10 


> 100 


> 20 kcV 


(McV) 





99.6 % 


97.0 % 


78.9 % 


91.6 % 


3.7 


0.001 


96.7 % 


71.9 % 


4.8 % 


29 % 


0.3 


0.01 


87.8 % 


59.5.6 % 


6.8 % 


22.8 % 


0.6 


0.1 


92.7 % 


50.4 % 


10.8 % 


20.9 % 


0.7 


0.5 


80.4 % 


60.4 % 


22.8 % 


39.6 % 


1.8 


1.0 


89.3 % 


72.2 % 


29.5 % 


44.5 % 


1.3 


5.0 


99.1 % 


92.8 % 


59.3 % 


78.6 % 


2.1 


10.0 


98.3 % 


93.2 % 


66.8 % 


85.0 % 


1.8 


20.0 


98.7 % 


95.8 % 


68.1 % 


89.6 % 


1.5 


50.0 


98.9 % 


97.2 % 


87.2 % 


94.9 % 


3.2 


100.0 


78.3 % 


68.4 % 


49.8 % 


65.3 % 


2.2 


500.0 


54.3 % 


50.6 % 


20.8 % 


43.9 % 


0.6 


1000.0 


41.3 % 


33.0 % 


3.9 % 


23.2 % 


0.3 


10000.0 


0.1 % 


0.00 % 


0.00 % 


0.00 % 


0.005 


1 < 1 1 f 1 V : ' ) . 




ye of Accelerated Electrons for Eq = 


U.UUUi 














> 1 


> 10 


> 100 


> 20keV 


(MeV) 





98.7 % 


88.5 % 


0.0 % 


50.4 % 


0.43 


0.001 


75.5 % 


31.6 % 


0.0 % 


3.4 % 


0.07 


0.01 


76.4 % 


21.7 % 


0.6 % 


2.8 % 


0.12 


0.1 


72.5 % 


22.1 % 


1.9 % 


2.0 % 


0.097 


0.5 


60.6 % 


31.9 % 


4.3 % 


7.5 % 


0.18 


1.0 


65.5 % 


34.7 % 


5.0 % 


9.3 % 


0.13 


10.0 


63.0 % 


32.4 % 


4.0 % 


6.9 % 


0.1 


20.0 


61.9 % 


38.4 % 


4.7 % 


7.7 % 


0.082 


50.0 


84.6 % 


62.0 % 


11.2 % 


23.0 % 


0.28 


100.0 


60.2 % 


35.0 % 


3.6 % 


7.2 % 


0.11 


500.0 


39.5 % 


11.7 % 


0.5 % 


0.1 % 


0.036 


1000.0 


30.3 % 


6.7 % 


0.2 % 


0.0 % 


0.022 



tion region: protons that can be accelerated at all no longer 
gain much energy after this time. 

4. Discussion and conclusions 

In this work we investigate the particle acceleration conse- 
quences of time-varying electric fields superposed on a X- 
type magnetic field to mimic generic features of dynamic, 
coUisionless reconnection (Craig and McClvmont 11991] ). We 
have shown that protons and electrons may gain relativis- 
tic energies in times <1 s; for plausible (small) electric field 
amplitudes and active region magnetic fields. This para- 
metric study is meant to be complementary to Petkaki and 
MacKinnon (|1997p where we attempted to tie this test- 
particle approach self-consistently to an MHD description 
of the passage of a wave. Although there are some quali- 
tative similarities with the particle behaviour studied ana- 
lytically by Litvinenko (|2003p . his adoption of a spatially 
uniform electric field makes direct comparison difficult. 

Before discussing some consequences of our results, we 
note some limitations of our calculation. First, this is a 
test particle approach. Particles do not interact with each 
other, nor do they influence the background field. In par- 
ticular, the particle distribution including the accelerated 
component may well be unstable to growth of various sorts 
of waves. Obviously such wave growth would influence the 
motion of particles, but we neglect this possibility. We ne- 
glect also radiation losses. In the solar corona this is not 
a serious neglect (even for 10 MeV electrons the radiative 
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energy loss time is 3000 s), but elsewhere in the cosmos 
it could become significant. 

We followed Craig and McClymont (|199ip in assuming 
a smoothly varying X-type field through the whole of our 
system. We took the active region lengthscale of 10® cm 
to define the boundary of the system. It is possible that 
the field strength increases more rapidly from the neutral 
point, approaching a constant value at smaller distance. 
This would reduce both the adiabaticity radii of particles 
and the characteristic timescale. If we continued to scale 
the electric field region with the proton adiabaticity radius, 
nothing would change except that the integration periods 
correspond to smaller real times. Thus particle acceleration 
would proceed more rapidly; however, fewer particles would 
be involved. Further consideration of this question might 
proceed via study of more realistic configurations includ- 
ing a neutral point or sheet (e.g. Forbes and Priest |1995( 
Fletcher & Martens HMEi Titov and Demoulin HMD or 
complex magnetic field structures (e.g. Malara et al. l2000| . 

The finite width of the nonadiabatic region allows par- 
ticles to gain or lose some energy randomly before return- 
ing to adiabatic motion. Together with repeated encounters 
with the dissipation region, the consequence of mirrorings 
in the extended configuration, this results in a Fermi-type, 
'stochastic' acceleration. 

In our model particle acceleration takes place for ge- 
ometrical reasons. The test particle calculation is numer- 
ically simpler than self-consistent approaches (e.g. Vlasov 
simulations, see Petkaki et al. (|2003l I2006P and gives use- 
ful insights to the particle energization process. There is no 
threshold for this type of acceleration, unlike resonant inter- 
action with low- frequency, MHD waves. The necessity for 
protons particularly to have threshold energies of around 
25 keV is a well known difficulty when such mechanisms 
are invoked (e.g. Forman et al., [1986). Our results indi- 
cate that low-frequency waves may themselves perform the 
'first-step' acceleration, if they propagate in a coronal struc- 
ture including a neutral point. This may occur indepen- 
dently of, or simultaneously with, the resonant cascade sce- 
nario of Miller and Vinas (|1993p . Possible difficulties with 
the number of pre-accelerated particles may be obviated if 
many neutral points are present, although such a situation 
obviously needs separate investigation (KliemT994|). 

Most of the resulting proton distributions have a bi- 
modal form (see Fig. [3]). Electron distributions are also 
bi-modal for the highest frequencies, 20 < w < 500 (see 
Fig. [S]). Whereas for the lowest frequencies of the electric 
field the bulk of the initial electron Maxwellian distribution 
is accelerated, for the highest frequencies only part of the 
electron distribution is accelerated (Table S]) . Acceleration 
occurs for all frequencies w < 10 when addressing the pro- 
ton distributions (Table [T]) . The bi-modal form of the pro- 
ton energy distributions might offer a way to have protons 
of gamma-ray producing energies {Kkin ~ 2 MeV) with- 
out the energetically dominant population at lower ener- 
gies that is the inevitable consequence of a diffusive par- 
ticle accelerator (see the Appendix of Eichler (1979); and 
MacKinnon (1991)). Investigation of the velocity space sta- 
bility of these distributions needs details of the angular dis- 
tribution at particular points in space, and is not discussed 
here. 

We note the effectiveness of acceleration of the two 
species varies according to the frequency of oscillation in- 
voked. Electrons are accelerated for a broader spectrum of 



frequencies. Frequencies 0.001 to 1000 have been simulated 
here, corresponding to real frequencies in the range 5 Hz to 
5 MHz (cf. the frequency range of waves from the base of 
the solar corona, probably in the range 0.01 Hz to lOKHz, 
e.g. Marsch et al. (|1982p '). Frequencies lower than 0.001 will 
also accelerate electrons as indicated by the net accelera- 
tion achieved for the constant electric field cases (Fig. [7|) , 
but frequencies higher than 1000 do not produce a net ac- 
celeration in the timescale of our model. 

Considered as a function of uj, the mean energy of the ac- 
celerated electron distribution exhibits a peak in the broad 
range 5 < w < 100. Such a peak leads us to suspect a reso- 
nance involving two or more of the timescales in the prob- 
lem. The initial gyrofrequencies of electrons lying in the 
adiabatic portion of the dissipation region also generally lie 
in this range. Inverse crossing times {1/tcr, see Eguation fH)) 
comparable with to might also lead to enhanced accelera- 
tion. Using Equation [Ml but taking account also of the 
mean increase in u^.y we do indeed find upper limits in the 
range 5 < 1/tcr < 100. 

Protons are accelerated for low electric field frequen- 
cies, achieving 7-ray producing energies in 5360 Tp = 1 s 
for frequencies w < 10 and for Eq = 0.001 and Eq = 0.01. 
A local peak in the mean energy of the accelerated proton 
distribution is seen at 0.1 < uj < 2.0. This range of fre- 
quencies are comparable to the gyrofrequencies of protons 
in the adiabatic region for our set of initial conditions and 
to the proton inverse crossing time. 

The variability of the effectiveness of acceleration of 
the two species according to the frequency of electric field 
oscillation might bear on the apparent variation of elec- 
tron/proton ratios in fiares (Ramaty & Murphy '1987') and 
the phenomenon of 'electron-only' flares (Rieger 1989). As 
a general comment, we note that higher frequency distur- 
bances favour electrons over ions, although more definitive 
statements will need a proper treatment involving a more 
realistic wave. 

For most frequencies and for constant electric field, part 
of the electron distribution escapes from the system bound- 
aries before 230400Te = 1 s. Electrons on average escape in 
less than 0.6 seconds in the frequency range 10 < a; < 50. 
Protons on the other hand do not escape the system bound- 
aries on the same timescale of 5360Tp = 1 s for Eq = 0.001 
and Eq = 0.0001. Electrons are accelerated more rapidly 
than protons to energies that do not allow them to mirror 
inside our system boundaries. 

Here and in Petkaki & MacKinnon (|1997p we investigate 
particle acceleration at a null in the presence of a linear 
disturbance. Such calculations may give some insight into 
particle acceleration in flares, although conditions then pre- 
sumably depart severely from linearity, but might be most 
relevant to quiescent, long-lasting phenomena such as ra- 
dio noise storms. Definitely involving deka-keV electrons 
(Raulin & Klein, '1994') and showing correlations with X-ray 
variations, but without chromospheric, fiare-like signatures 
(Svestka et al 119861 Crosby et al. [TMBj), particle accelera- 
tion in noise storms might occur as described here, if the 
relevant coronal structures include null points. Electrons 
accelerated at a neutral point will likely encounter very 
large mirror ratios, trapping them in the corona (Fletcher 
&; Martens [1998]) and accounting for the exclusively coronal 
phenomena accompanying noise storms. 
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Fig. 4. Mean Energy of proton distributions for three am- 
plitudes of the electric field {Eq ^ 0.0001,0.001,0.01) and 
for range of frequencies 0.001 to 100.0. The mean energy 
of the initial Maxwellian distribution is shown as a straight 
full line. The constant electric field case is represented by 
uj = 0.0001. 
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Fig. 5. Standard deviation of proton distributions for three 
amplitudes of the electric field {Eq = 0.01,0.001,0.0001) 
and for range of frequencies 0.001 to 100.0. The standard 
deviation of the initial Maxwellian distribution is shown 
as a straight full line. The constant electric field case is 
represented hy uj — 0.0001. 



Fig. 2. Proton orbit in time- varying electric field of w = 0.2. 
(a) Energy as a function of time (b) Projection in the X-Y 
plane, (c) Projection in the X-Z plane. 
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Fig. 7. Mean Energy and standard deviation of electron 
distributions for two amplitudes of the electric field (Eq = 
0.001,0.0001) and for range of frequencies. The mean en- 
ergy of the initial Maxwellian distribution is shown as a 
straight full line. The constant electric field case is repre- 
sented in a; = 0.001. 




Fig. 8. Time evolution of logarithm of mean energy of pro- 
tons for amplitude of the electric field {Eq = 0.0001) for 
constant electric field (represented on the graph by 0.0001) 
and for range of frequencies 0.001 to 100.0. 



